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We present uniaxial pressure structural relaxations for CaFe2As2 and BaFe2As2 within density 
functional theory and compare them with calculations under hydrostatic pressure as well as available 
experimental results. We find that CaFe2As2 shows a unique phase transition from a magnetic 
orthorhombic phase to a nonmagnetic collapsed tetragonal phase for both pressure conditions and 
no indication of a tetragonal phase is observed at intermediate uniaxial pressures. In contrast, 
BaFe2As2 shows for both pressure conditions two phase transitions from a magnetic orthorhombic 
to a collapsed tetragonal phase through an intermediate nonmagnetic tetragonal phase. We find 
that the critical transition pressures under uniaxial conditions are much lower than those under 
hydrostatic conditions manifesting the high sensitivity of the systems to uniaxial stress. We discuss 
the origin of this sensitivity and its relation to experimental observations. 

PACS numbers: 74.70.Xa,61.50.Ks,71.15.Mb,64.70.K- 



The discovery of superconductivity in iron pnictides [T] 
has initiated an enourmous amount of activities related 
to these materials. Superconductivity can be triggered 
either by chemical doping or by application of pressure 
on the undoped parent compounds. One of the fam- 
ilies that has been intensively studied under pressure 
is the 122 family AEFe 2 As 2 (AE = Ca, Sr, and Ba). 
CaFe2As2 at ambient pressure undergoes a first order 
phase transition from a tetragonal to an orthorhombic 
phase at 172 K accompanied by a magnetic transition. 
Initial reports on pressure experiments showed that at 
P ~ 0.23 GPa the orthorhombic and antifcrromagnetic 
phases are suppressed and the system superconducts at 
low temperatures [3] [3]. Moreover, a compressed tetrago- 
nal phase - also called 'collapsed' tetragonal phase - was 
identified at higher pressures. Subsequent susceptibil- 
ity and transport measurements under hydrostatic con- 
ditions showed at low temperatures and P ~ 0.35 GPa 
a sharp orthorhombic to collapsed tetragonal phase but 
no signature of superconductivity [4 . In contrast, recent 
neutron diffraction experiments on CaFe2As2 under uni- 
axial pressure along the c axis [5] indicate for pressures 
above 0.06 GPa and low temperatures the presence of an 
intermediate nonmagnetic tetragonal phase between the 
magnetic orthorhombic and the nonmagnetic collapsed 
tetragonal phases. This phase was identified by the au- 
thors as the phase responsible for superconductivity at 
T = 10 K. Other reports based on muon spin-relaxation 
measurements suggest the existence of superconductivity 
in the orthorhombic phase, raising the question whether 
superconductivity and magnetism can coexist [§]. 

BaFc2As2 shows an even more complex behavior under 
pressure. At ambient pressure it undergoes a phase tran- 
sition from a metallic tetragonal phase to an orthorhom- 
bic antifcrromagnetic phase at T = 140 K. Under pres- 
sure the gradual appearance of a superconducting dome 



has been observed by various groups [7J though the role of 
nonhydrostatic conditions is not yet well understood [S]. 
Recent synchrotron X-ray diffraction experiments under 
pressure [3] observe at T = 300 K a tetragonal to col- 
lapsed tetragonal phase transition at P — 22 GPa under 
hydrostatic conditions while this transition appears al- 
ready at P = 17 GPa under nonhydrostatic conditions. 
On the other hand, the authors of Ref. [TO] find at a lower 
temperature of T = 33 K that BaFe2As2 undergoes a 
phase transition from a magnetic orthorhombic to a non- 
magnetic collapsed tetragonal phase at P = 29 GPa and 
report an anomaly in the As-Fe-As bond angles at 10 GPa 
that they ascribe to be of electronic origin. In contrast, 
high-pressure neutron diffraction experiments [llj per- 
formed at T — 17 K find a tetragonal phase at 3 GPa 
and 6 GPa. 

These experimental results show that the onset of 
superconductivity as well as the appearance of several 
structural phases at low temperatures in CaFe2As2 and 
BaFe 2 As 2 are extremely sensitive to the pressure condi- 
tions [5J [T2Trr4] and are a subject of intensive debate. In 
view of this strong controversy we present ab initio den- 
sity functional theory results for the electronic, magnetic 
and structural behavior of both systems under uniaxial 
and hydrostatic pressure conditions. Previous theoreti- 
cal approaches which have examined the properties of the 
122 family under hydrostatic pressure have employed ei- 
ther fixed volume structural optimizations [15] or molecu- 
lar dynamics [16) . Recently, anisotropic pressure studies 
on BaFe 2 As 2 based on ground state geometry calcula- 
tions of more than 300 structures at different fixed vol- 
umes were reported in Ref. [17] . Our approach consists 
of constant pressure structural relaxations which allows 
us to probe the low-temperature portion of the phase 
diagram in a relatively simple and straightforward way. 
With this approach we can treat nonhydrostatic condi- 
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tions which are at the heart of this work. 

We find that uniaxial pressure along the c axis sig- 
nificantly reduces the transition pressures in both sys- 
tems. CaFe2As2 shows for both pressure conditions an 
orthorhombic to collapsed tetragonal transition, though 
the transition is less abrupt when uniaxial pressure is 
applied. For BaFe2As2 we observe two phase transitions 
from orthorhombic to collapsed tetragonal through an 
intermediate nonmagnetic tetragonal phase. An analysis 
of the electronic bandstructure features near the criti- 
cal pressures reveals the origin of the sensitivity of the 
systems to pressure conditions. 

Calculations were performed using the Vienna ab ini- 
tio simulations package (VASP) [18] with the projector 
augmented wave (PAW) basis [19] in the generalized gra- 
dient approximation (GGA). Structural relaxations un- 
der hydrostatic pressure were carried out with the conju- 
gate gradient (CG) method as implemented in the VASP 
package. The energy cutoff was set to 300 eV and a 
Monkhorst-Pack uniform grid of (6 x 6 x 6) points was 
used for the integration of the irreducible Brillouin zone. 
For relaxations with the CG algorithm two cycles were 
performed in order to minimize the error caused by the 
Pulay stress. Note that the reported bond compressions 
of up to 7% at 50 GPa don't affect the precision of the 
PAW basis. In order to perform relaxations under uni- 
axial pressure we modified the fast inertial relaxation 
engine [20] (FIRE) algorithm to be able to handle full 
structural relaxations with an arbitrary stress tensor. 

In Fig. [l]we show the evolution of lattice parameters, 
volume and Fe-Fe, Fe-As distances under hydrostatic and 
(c-axis) uniaxial pressure for CaFe2As2- We find a first 
order phase transition from a magnetic (stripe order) or- 
thorhombic phase to a nonmagnetic collapsed tetragonal 
phase at P c = 3.05 GPa (P c = 0.48 GPa) under hy- 
drostatic (uniaxial) pressure and at zero temperature we 
don't observe any intermediate tetragonal phase under 
uniaxial stress [5]. 

In hydrostatic conditions, a and b expand as a conse- 
quence of the Pauli principle 16 with b abruptly increas- 
ing in value, while c shows a significant collapse of 6.5% 
(Fig. [I] (a)) and the unit cell volume shows a sharp drop 
of about 4.3% (Fig.[T](c)). The value of c/a t = 2.58 with 
at = a/v2, indicates the onset of a collapsed tetrago- 
nal phase. Our results are in good qualitative agreement 
with experimental [3] observations, except for the over- 
estimation of the critical pressure (-P,? xp = 0-3 GPa) also 
found in previous theoretical studies [HI [16] . Following 
the changes of the lattice parameters at P c , the inplane 
Fe-Fe distances show a sharp increase at P c while the out- 
of-plane Fe-As distance decreases (Fig. [I] (e)). Using the 
generalized Birch-Murnaghan p—V equation of state [21] 
we obtained a bulk modulus B — 70 ± 3 GPa at ambient 
pressure, while at P c the bulk modulus jumps from 56 ±3 
to 105 ± 2 GPa. In order to obtain these estimates we 
performed a series of fits for every phase separately con- 



sidering every pressure point of our data as a reference 
pressure. In this way we obtain the bulk modulus as a 
function of pressure. 

In contrast to the hydrostatic case, when uniaxial pres- 
sure is applied (Fig. [T] (b)) the a and b lattice param- 
eters expand significantly while c is compressed up to 
P c = 0.48 GPa where a drop for c is observed while a 
and b continue to expand monotonously. The volume re- 
duces by 3.4% and the ratio c/a t = 2.56 at P c (Fig.[l](c)) 
denotes the entrance to a collapsed tetragonal phase, 
where magnetism is suppressed completely. The phase 
transition shifts to smaller P c compared to hydrostatic 
pressure, which is in very good agreement with experi- 
ments under nonhydrostatic pressure conditions [5 ; . Nev- 
ertheless, the authors of Ref. [5] find for pressures above 
0.06 GPa a stabilization of the high-temperature tetrago- 
nal structure down to temperatures below the supercon- 
ducting transition. This phase is not seen in our calcula- 
tions which may be related to the fact that at very low 
temperatures the tetragonal phase may be disappearing 
again (Fig. 1 of Ref. [5]). 

In order to understand the differences in behavior ob- 
served between the hydrostatic and uniaxial pressures, 
we show in Fig. [2] the orbital weighted bandstructure 
and k z = 0, k v = and k x + k y =0 Fermi surface 
cuts of CaFe2As2 under hydrostatic (Fig. [2] (b)-(c)) and 
uniaxial pressure (Fig. [2] (d)-(e)) at pressures below and 
above the phase transition. Bandstructures and Fermi 
surfaces were calculated using the full-potential local or- 
bital (FPLO) basis 023- The bandstructure and Fermi 
surface cuts at ambient pressure are also shown for com- 
parison (Fig. [2] (a)). We use the orthorhombic space 
group F mmm for all band structure plots in order to fa- 
cilitate comparison. The behavior of the electronic struc- 
ture in the vicinity of the Fermi energy is crucial for un- 
derstanding the transition. Right below P c both pressure 
conditions show a high density of Fe d xz , d yz and d x 2_ y 2 
states at E F (see in Fig. [2] (b) and (d) the T-Z path and 
near M) which are pushed away above P c (Fig. [2] (c) and 
(e)) and the hole pockets at V disappear, suppressing 
possible nesting conditions. The compression along c en- 
forces the interlayer As p z -As p z bonding [23j which can 
be related to the proximity of the As p z band to Ep near 
P c . Uniaxial stress is for this process more effective than 
hydrostatic pressure since a similar electronic behavior is 
reached at much smaller pressures as observed in Fig. [2] 
(d)-(e). In agreement with Ref. [24] the collapsed tetrago- 
nal phase sets in as soon as the Fe magnetic moment goes 
to zero. Also, the shape of the Fermi surface derived from 
Fig. [2] in the collapsed tetragonal phase agrees well with 
the de Haas van Alphen measurements performed for 
CaFe2P2 (c/a t — 2.59) where a highly dispersive topol- 
ogy in the c axis as well as the absence of the hole pocket 
at the T point has been reported [25] (compare Fig. [2] (c) 
and [2] (e) with Figs. 2 and 3 of Ref. [25]). The isoelec- 
tronic substitution of As by P in CaFe2As2 corresponds 
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FIG. 1. Structure of CaFe2As2 under hydrostatic and uniaxial pressure. Lattice parameters, a) and b), volume and axis ratio, 
c) and d), c) and d), selected bond lengths e) and f). 



to application of chemical pressure and shows similar fea- 
tures to the collapsed tetragonal phase of CaFe2As 2 ob- 
tained after application of (hydrostatic) pressure. The 
similarity of chemical pressure and applied pressure has 
already been discussed in Refs. f2"5l and ITT1 In fact, com- 
parison of our obtained c/at = 2.58 ratio and As position 
zas = 0.1358 (hydrostatic) in the collapsed tetragonal 
phase of CaFe2As2 with the measured c/a t — 2.59 and P 
position zp = 0.1357 of CaFe2P2 shows the high resem- 
blance between both crystal structures. 

We now proceed with the analysis of BaFe2As2. In 
Fig. [3] we present the changes in lattice parameters, vol- 
ume and atomic distances under hydrostatic and uniaxial 
pressures for BaFe2As2- Similar to CaFe2As 2 , the critical 
pressures under uniaxial stress are reduced with respect 
to hydrostatic conditions. This observation was also re- 
ported by recent constant volume density functional the- 
ory calculations on BaFe 2 As2 under nonhydrostatic pres- 
sure [T7]. BaFe2As2, contrary to CaFe2As2, shows two 
phase transitions. At P C1 = 11.75 GPa (P Cl = 0.72 GPa) 
we find a phase transition from an antiferromagnetic 
orthorhombic to a nonmagnetic tetragonal phase un- 
der hydrostatic (uniaxial) conditions. A second smooth 
phase transition to a collapsed tetragonal phase is ob- 
tained for P C2 = 28.6 GPa (P C2 = 3.17 GPa) (Fig. § 
(c)-(d)) [9J. High-pressure neutron diffraction experi- 
ments jllj as well as previous theoretical calculations un- 
der hydrostatic pressure conditions also find a phase tran- 
sition to an intermediate tetragonal phase [E3 [TB] but re- 
cent syncrotron X-ray diffraction experiments under non- 
hydrostatic conditions see no signature of an intermediate 



tetragonal phase at low temperatures. Nevertheless, an 
anomaly in the As-Fe-As bond angles at P ~ 10 GPa [TU] 
as well as a loss of magnetic moment [Uj have been re- 
ported. This could be related to the phase transition 
that we find at P C1 = 11.75 GPa where magnetism is sup- 
pressed. At higher pressures the agreement of the onset 
of the collapsed tetragonal phase with the X-ray diffrac- 
tion data [TU] is very good. Clearly the phase transitions 
in BaFe2As2 are less abrupt than in CaFe2As2- The am- 
bient pressure bulk modulus is estimated at 67±4 GPa, in 
good agreement with experimentally reported values |10j 
of 82.9±1.4 and 65.7±0.8 GPa at 33 K and 300 K, respec- 
tively. At P Cl , the bulk modulus abruptly increases from 
98 ±4 to 128 ±3 GPa and at P C2 it jumps from 150 ±3 to 
173 ±2 GPa. This is in very good agreement with the ex- 
perimental estimate of B = 153 ± 3GPa for the collapsed 
tetragonal phase [10]. We also analyzed the Fe-As bond 
compressibility for P = 9 GPa (hydrostatic) and found 
k = 3.5 x 10~ 3 GPa -1 which is in excellent agreement 
with k = 3.3 x 10~ 3 GPa -1 obtained in an extended X- 
ray absorption fine structure (EXAFS) experiment |26j . 
In Fig. g] we show the comparison of the measured pres- 
sure dependence of the Fe-As bond distances [26] with 
our results. Due to different temperatures (experiment 
is performed at room temperature, theory at T — 0) our 
distances are shorter by about 0.02 A (0.8%), but the 
overall agreement is good. 

In Fig. [5] we present the orbital weighted bandstruc- 
ture and k z — 0, k y — and k x + k y = Fermi surface 
cuts of BaFe2As 2 under hydrostatic (Fig. [5] (b)-(c)) and 
uniaxial pressure conditions (Fig. [5] (d)-(e)) at two pres- 
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FIG. 2. Bandstructure and k z = 0, k y = and k x + k y = Fermi surface cuts of CaFe2As2 [22) . For the orbital character, x 
and y point along the nearest neighbour Fe-Fe connections. 



sures below and above the orthorhombic to tetragonal 
phase transition at P Cl = 11.75 GPa (P Cl = 0.72 GPa). 
Similar to CaFe2As2 we observe below P Cl a high density 
of Fe d xz , d yz and d x 2_ y 2 states at Ep which is pushed 
down (less drastically than in CaFe2As2) for pressures 
above P C1 . The hole pockets disappear at the T point 
and the Fe magnetic moment goes to zero. Here the As 
p z band seems to be little affected at the critical pres- 
sure. In contrast, at P C2 = 28.6 GPa {P C2 = 3.17 GPa) 
(bandstructure not shown) the As p z band is pushed to- 
wards the Fermi level indicating a strong As p z -As p z 
bonding while entering the collapsed tetragonal phase. 
These results show that under perfect hydrostatic or per- 



fect uniaxial pressure conditions neither the intermediate 
tetragonal phase nor the collapsed tetragonal phase ful- 
fill Fermi surface nesting conditions. In fact, we find that 
the structural parameters measured in Ref. [TT] are sim- 
ilar to our calculated parameters far below P Cl in the 
orthorhombic phase (except for the orthorhombic distor- 
sion), where well defined hole pockets are found at the T 
point. 

In summary, we have presented finite pressure den- 
sity functional theory calculations which allow the in- 
vestigation of nonhydro static pressure conditions. Our 
finite pressure relaxed structures show good agreement 
with the available experimental data (volume, Birch- 
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FIG. 3. Structure of BaFe2As2 under hydrostatic and uniaxial pressure. Lattice parameters, a) and b), volume and axis ratio, 
c) and d), selected bond lengths e) and f). 
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FIG. 4. Comparison of T = calculated (squares) and T = 
298 K measured (triangles, circles) [26] Fe-As bond distances 
in BaFe2As2 as a function of pressure. 



Murnaghan values and compressibilities) though our 
magnetic moments in the orthorhombic phase are larger 
than the observed experimental values. We expect that 
this overestimation affects mostly the values of P c . Com- 
parison of our calculated Fe-As bond distances at differ- 
ent pressures with measured distances [26] shows good 
agreement. Also, available de Haas van Alphen measure- 
ments performed for CaFe2P2 [25] agree well with our 
predicted Fermi surface shapes for CaFe2As2 in the col- 



lapsed tetragonal phase. This overall agreement with ex- 
perimental observations demonstrates that the presented 
constant pressure calculations provide a reliable theoret- 
ical prediction of structures under nonhydrostatic pres- 
sure conditions, allowing for arbitrary stress tensors in 
future studies. Such calculations can complement ex- 
periments and help identify the precise degree of hydro- 
staticity. We find that uniaxial stress along the c axis 
considerably reduces the critical pressures for CaFe2As2 
and BaFe2As2- This behavior can be understood by the 
fact that the phase transitions are strongly dictated by 
the electronic properties in the vicinity of the Fermi en- 
ergy, as shown by our electronic structure analysis. While 
CaFe2As2 undergoes a magnetic orthorhombic to a non- 
magnetic collapsed tetragonal phase for both pressure 
conditions and no indication of an intermediate tetrag- 
onal phase is observed under uniaxial stress, BaFe2As2 
shows two phase transitions from a magnetic orthorhom- 
bic to a collapsed tetragonal phase through an interme- 
diate nonmagnetic tetragonal phase for both pressure 
conditions. All nonmagnetic phases show a disappear- 
ance of the hole pockets at the T point suppressing pos- 
sible Cooper pair scattering channels between electron 
and hole pockets. Such scattering channels have been 
discussed to be important for the superconductivity in 
BaFe2As2 [57]. More experiments need to be done in 
order to understand the origin of the superconducting 
phase observed in these materials under various pressure 
conditions. 
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